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Abstract 

Generalizing recent work, the Raman scattering intensity from a semi-infinite 
superconducting superlattice is calculated taking into account the surface 
contribution to the density response functions. Our work makes use of the 
formalism of Jain and Allen developed for normal superlattices. The sur- 
face contributions are shown to strongly modify the bulk contribution to the 
Raman-spectrum line shape below 2A, and also may give rise to additional 
surface plasmon modes above 2A. The interplay between the bulk and sur- 
face contribution is strongly dependent on the momentum transfer qn parallel 
to layers. However, we argue that the scattering cross-section for the out-of- 
phase phase modes (which arise from interlayer Cooper pair tunneling) will 
not be affected and thus should be the only structure exhibited in the Raman 
spectrum below 2A for relatively large qn ~ OTA/up. The intensity is small 
but perhaps observable. 

PACS numbers: 74.20.Fg, 71.45. Gm, 74.80.Dm, 78.30.-j 
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I. INTRODUCTION 



Recently, the authors [T|] studied the inelastic light-scattering intensity of a semi-infinite 
superconducting superlattice with a bilayer basis. Motivated by the Cooper pair tunneling 
model proposed by Chakravarty, Anderson, and coworkers for high T c layered super- 
conductors, we discussed the in-phase and out-of-phase phase modes (corresponding to the 
phase fluctuations of the two superconducting order parameters in a bilayer) which arise in 
the presence of interlayer Cooper pair tunneling ||. These modes couple into density fluc- 
tuations and, as a result, show up in the Raman inelastic light scattering. The intensity is 
weak because of screening associated with the Coulomb interaction, and in fact is below the 
threshold of current Raman experiments. However, our results are of sufficient interest that 
such experiments should be attempted. In this paper, we extend our previous calculations 
HJ and show how the Raman spectrum for uj < 2A is significantly modified when we include 
the surface contribution. It is also the first time surface plasmons (which occurs above 2A) 
are included in the calculation of the Raman intensity of layered superconductors. 

The present paper is based on the approach of Jain and Allen ||, who considered nor- 
mal layered electron gas (LEG). In their calculation for a semi-infinite superlattice, both 
the bulk and surface contributions were included. They found that there were two effects of 
the surface: (1) Van Hove singularities at the upper and lower limits of the bulk plasmon 
band were completely canceled out by negative surface contributions. (2) Depending on the 
background dielectric constants, surface plasmons || can appear, either above or below the 
bulk plasmon band. We show that the analogous effects arise in a semi-infinite supercon- 
ducting superlattice, resulting in major modifications of the bulk contribution to the Raman 
scattering spectrum given in Ref . |l| . Apart from the out-of-phase phase mode contribution, 
the Raman scattering intensity is found to be strongly dependent on the value of momentum 
transfer q\\ (parallel to the layers). For simplicity, we only discuss the surface effects for a 
superlattice with a single layer per unit cell. This is sufficient to understand the essential 
physics. 
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The isotropic inelastic light-scattering cross section is given by |4],^] 

d(T A i2 7/ \ /-,\ 

— Qc|e i -e / | I{q, U ), (1) 

where 

I(q,u,) =^lmXiMbuJ,l'>~ iZt>i+Zl ' j)/S 
l,l' 

X e -*k±(.Zl,i-Z t >j), (2) 

Here represents the position of the z-th layer in the l-ih unit cell. The density response 
function Xiji^b u i U i n (i) nas Deen evaluated in Ref 0] and represents the correlation 
between the charge density on layer and the charge density on layer The inci- 

dent photon has momentum kj, energy u>i, and polarization e; and the scattered photon is 
similarly described by kj, Uf, and e/. We assume that the energy transfer uo = uo,i — uj f in 
(0) is very small compared to the photon frequencies, i.e., cjj ~ Uf. The momentum transfer 
parallel to the interface is kiu — kn j = qii and the momentum perpendicular to the interface 
is k z> i — k z j = q z . For small-angle scattering, we have Reg z ~ 2k± and Img z ~ <5 _1 where 
k± is the momentum carried by the incident photon and 5 describes the damping of the 
photons in the medium. The result given in fl2|) shows that as a result of the finite value 
of S, the inelastic light-scattering cross section involves a weighted sum of the correlation 
functions for electronic densities in the different layers. We are interested in the interplay 
between the bulk and surface contributions to the Raman spectra. For this purpose, we only 
keep the isotropic matrix element for the Raman interaction given in ([!]) 0. We consider 
superconductors with s-wave layer pairing interaction but similar calculations could be done 
for e?-wave superconductors, as discussed in Refs. 0|§. 

II. BULK IN-PHASE AND OUT-OF-PHASE PHASE MODES 

For later comparison, we first recall the results of Ref. |1[ for the Raman spectra for a 
semi-infinite superconducting bilayer superlattice, ignoring the surface contributions. One 
finds that @ reduces to 
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where we have introduced the functions 



b = cosh q\\c — 2v2T>E + sinh q\\c 



u = b + Vb 2 -l 

F = u 2 e 2c/s - 2ue c/s cos 2k ± c + 1 
D_ = 1 - vWl - e^^L. 



(4) 



The spacing of the bilayer is d, the unit cell length is c, t>2o = 2vre 2 /g||e is the 2D Coulomb 
interaction, and e is the superlattice background static dielectric constant. In the long- 
wavelength limit (q\\ <C 2A/vp), the functions E±(q\\,u) are given by 



E± 

where we have defined 
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gN(e F ) x 2 - 1 

Here oj = lu/2A, q» = q\\v F /2A, N(ep) = m* /n is the 2D electronic density of states at the 
Fermi level with m* being the effective electronic mass, g is the in-layer pairing interaction, 
and Tj is the interlayer Cooper pair tunneling strength. Replacing oj — > oj + ry is a simple 
way of including finite energy-resolution. On the rhs of (|3|), the first term (= I/) gives the 
contribution from the in-phase phase fluctuations, while the second term (= Io) is associated 
with the out-of-phase phase fluctuations. 



One finds that the in-phase first term in (|3|) has three poles, given by 

F(2k±,q^,u)=0 and b(q h uj) = ±l. (8) 

F = gives an in-phase plasmon mode which Raman scattering picks up (in an approximate 
way 0], this mode is similar to the plasmon mode of an infinite superlattice, with q z = 2k±). 
The additional two (Van Hove) singularities given by the solutions of b = ±1 correspond to 
the upper (+) and lower (— ) limits of the "bulk plasmon" band for an infinite superlattice 
H]. In contrast, the second term in (||) only has a single pole given by 

D_(q h u)=0, (9) 

corresponding to out-of-phase phase mode discussed in detail in Ref. [|]]. Because the unit 
cell summation is over many bilayers, the Raman intensity is strongly enhanced in the 
superlattice case compared to the isolated bilayer case (the latter is discussed in Ref. ||). 
This is the origin of the prefactor [1 — exp(— 2c/8)]~ 1 ~ 5/2c in (^). Using c = 12 A and 
5 = 1000A, this prefactor is ~ 40. 

We note in (|3|) that, in the limit of q\\ — > 0, we have R + = and hence E + ~ g| — > 0. This 
implies that the in-phase phase modes given by F = and b = ±1 have less weight when 
qu is small, an expected consequence of the screening due to the Coulomb interaction. In 
contrast, even in the low-^n limit, E_ ~ R_ ~ x is still finite, being proportional to the pair 
tunneling strength Tj. This means that the out-of-phase phase mode given by D_ = has 
a weight proportional to x and is not too dependent on the value of q\\ (in the range probed 
in Raman scattering experiments). In addition, one can see from (|3|) that the intensities 
of // and Io are also dependent on the factors (1 + e~ 2d / s ± 2e~ d ^ 5 cos 2k±d) which arise 
from the lattice summation in (Q). Since d 5 and k^d < 1, we see that the intensity for 
the out-of-phase phase modes (— sign) is greatly reduced compared to that of the in-phase 
phase modes (+ sign) as a result of these factors. 

In Fig. [I], we plot the Raman light-scattering intensity based on ([3]). In this and other 
figures, we use the parameters: bilayer spacing d = 3A, the unit cell size c = 12A, pairing 
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strength gN(ep) = 0.25, Fermi momentum kp = 3.07 x 10 7 cm _1 and hence a 2D hole 
density n = 1.5 x 10 14 cm -2 , layer effective mass m* = m, background static dielectric 
constant e = 10, photon momentum in the ^-direction k± = 1.0 x 10 5 cm _1 , the optical 
penetration depth 5 ~ l/k± = 1000A, the superconducting energy gap A = 280cm _1 , 
and the finite-energy resolution 7 = 0.05A. The momentum transfer parallel to the layers 
is q\\ = 5.0 x 10~ 3 A/vf = 1.17 x 10 3 cm -1 in Fig. |l|. The out-of-phase phase mode is 
well-defined, as expected. While the Raman intensity from the out-of-phase contribution 
shown in Fig. [I] is not too dependent on the value of q\\, roughly speaking, the intensity 
from the in-phase contribution is proportional to gjj. As shown in Fig. [l] for the in-phase 
contributions, one has a small peak at u = 2A corresponding to the pair-breaking gap in 
an s-wave superconductor (which is identical to the upper limit of the bulk plasmon band, 
i.e., the pole given by b = 1). In addition, "hidden" in the low-frequency broadened peak 
is the in-phase phase mode contribution given by the solution of F = 0, which overlaps on 
the Van Hove singularity corresponding to the pole b = —1 at the lower limit of the bulk 
superlattice plasmon band (see Fig. |2|). This will become more transparent when we discuss 
the surface contribution. 

In a normal metal superlattice ||, the lower limit of the bulk plasmon band is far 
away from the particle- hole continuum and, as a result, both the upper (6 = +1 pole) 
and lower (b = — 1 pole) limits of the bulk plasmon band are well-defined. We recall that 
the bulk plasmon band refers to the plasmons labeled by q z in an infinite superconducting 
superlattice. These give rise to the Van Hove singularities discussed in Ref. [|J. In contrast, 
in a superconducting superlattice, the particle-hole excitation spectrum (which begins at 
the pair-breaking gap 2A) is strongly coupled into the superlattice bulk plasmon spectrum. 
As a consequence, the bulk plasmon band is split into two different regions above and below 
the pair-breaking gap (2A). For the plasmon band below 2 A (which we are most interested 
in), one can find a well-defined line for b = — 1 corresponding to the lower band limit of 
the bulk plasmon band (which is generally at low frequencies). However, due to the strong 
coupling between the bulk plasmon and BCS particle-hole continuum, starting at 2A, there 
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is no well-defined solution for b = +1. Nevertheless, the peak at 2 A in Fig. [T] supports the 
argument that u = 2A can be considered as the effective upper limit of a superconducting 



bulk plasmon band P,[T0||. 

As Jain and Allen [|J have pointed out for a normal superlattice with one layer per unit 
cell, the Van Hove singularity associated with b = +1 corresponds to all the neighboring 
layers oscillating in-phase; while the one associated with b = — 1 corresponds to all the 
neighboring layers oscillating out-of-phase with each other. In contrast, the out-of-phase 
phase mode in Fig. [l| is a collective mode associated with the "internal Cooper pair dynamics" 
exhibited by a bilayer via the interlayer Cooper pair tunneling between the two layers. The 
physics of this out-of-phase phase mode is completely different from the out-of-phase b = — 1 
bulk plasmon in a superlattice. 



III. SURFACE CONTRIBUTIONS AND SURFACE PLASMONS 

The result in does not include the surface contribution. We now include it but only 
consider the case of a semi-infinite superlattice of single layer per unit cell since this already 
describes the interplay between bulk and surface contributions. In the case of a single layer 
per unit cell, we need only to replace the usual 2D Lindhard function in the formulas given 
by Jain and Allen [|J by the appropriate density response function, [i.e., E + in (^)] for a 
neutral ZD superconductor 0. Using Eq. (50) in Ref. [Q, we find that the resulting Raman 
intensity is given by 

v 2dE + sinh q\\c(u 2 e 2c l s — 1) \ 



i + 



+ 



v 2 DE + {e 2c / s - l)(u 2 A - 2uB + C) 



2Q(b 2 - l)F 
where we have defined 



(10) 



A = Gsinh 2 g||C+l + |e 29 n c , 

B = if sinh 2 q\\c + coshgyc + 7^ el3 " C ! 

C = Gsinh 2 g|| C +l + |, (11) 



with 

G = -^[{b 2 — I) -5 — 1/ sinhgyc]/ sinhgyc, 

H = -[u~ x (b 2 — — e~ 9||C / sinhgyc]/ sinhgyc, 

Q = — [1 — (6 2 — — 6coshg||c)/ sinhgyc] 

— i«e 9 n c (6 2 — l)~5(coshg||C — 6)/ sinhgyc, (12) 

where w, and F are defined in (Q) and now b = coshgyc — t>2D-E+ sinhgiic. The parameter 
a = (e — e )/(e + e ) depends on the optical dielectric constants (e) inside and (e ) outside 
the superlattice. It plays a key role in determining the surface contributions as well as the 
appearance and energy of surface plasmons. The formula for a can be rewritten in the useful 
form e/eo = (l + a)/(l — a). We call attention to the similarity between the first term in (|10D 
and the first term in (|3|). In the rhs of (0), the first term (= i#) gives the bulk contribution, 
while the second term (= Ig) is associated with the surface contribution. The three poles 
mentioned earlier are exhibited by both contributions: F = corresponding to an in-phase 
plasmon which Raman scattering picks up and b = ± 1 corresponding to the upper and lower 
limits of the bulk plasmon band. There is a new pole of the surface contribution Ig given 
by 

Q(q h u)=0, (13) 

which can be shown to correspond to a surface plasmon. One has a nontrivial solution of 
Q = only when a^O (i.e., e ^ eo), that is, the surface of the superconducting superlattice 
must separate regions with different dielectric constants to give rise to surface plasmons. 

In Fig. |2|, we show the dispersion relation of surface plasmons in a superconducting 
superlattice, as given by the solutions of Q = 0, for various values of a. The shaded area 
between the line denoted by b = —1 to the line ui = 2 A represents the bulk plasmon band 
of an infinite superconducting superlattice. We find that the surface plasmon appears only 
above the upper limit (i.e., uo = 2A) of this bulk plasmon band, where BCS particle-hole 
damping of collective modes can occur. For positive values of a (e > eo), a surface plasmon 
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appears at very large energies ui ^> 2 A, in which case it is essentially identical to that in 
a normal superlattice. We remark that for a high-T c material with a dielectric constant 
e ~ 10 in a vacuum (eo = 1), one has a = 0.82. For a increasingly negative (— ► —1.0), which 
requires eo/e > 1, the surface plasmon energy slowly decreases toward 2 A (see Fig. |2|). 

In Fig. |2], the dispersion relation denoted by F = represents a (bulk) plasmon, which is 
a pole of the Raman scattering intensity in (|T0|) . One sees that F = mode is well-defined 
only when q\\ < 0.025A/-uf (solid line). The critical value of q\\ (0.025A/Vp) changes for 
different choices of 5 and k±. When q\\ > 0.025A/Vf, we find no solution for F = 0. The 
dashed line represents the minima of F (i.e., an over-damped or relaxational mode). This 
broad resonance is always peaked at u/2A = 0.8, whatever the values chosen for 5 and kj_. 

In Fig. [|, we show the net Raman intensity based on (|I0|) for a semi-infinite supercon- 
ducting single-layer superlattice, showing the surface and the bulk components. Comparing 
Fig. |^ with Fig. p], one sees that the low- weight Van Hove singularity at b = — 1 from the 
bulk contribution is canceled by the negative surface contribution. In contrast, the in-phase 
phase bulk plasmon mode (F = 0) shows up as a sharp peak in the low-frequency region. 
As shown in Fig. [| a peak associated with the 6 = 1 Van Hove singularity can still appear 
at u — 2A since the bulk and surface contributions do not completely cancel each other. 

Fig. £| is similar to Fig. [| using the same parameters, but at a much higher momentum 
transfer gn = O.lA/vp. The cancellation between the bulk and surface contributions is 
clearly shown not only for the two boundaries of the bulk plasmon band at b = — 1 and 
uj = 2 A, but almost for the entire region below 2A. We also see that at q» = O.lA/vp, the 
broad relaxational mode [always peaked at u « 0.8(2A)] corresponding to the minimum of 
F (see Fig. 0) has only low weight in the Raman scattering spectrum. 

To see how the interplay between the bulk and surface contribution depends on gn, we 
plot in Fig. |^ the total Raman intensities for various values of q» . Comparing Fig. |5] with 
Fig. 0, one finds that for q\\ % 0.025A/Vf, a well-defined bulk plasmon mode (F = 0) gives 
rise to a sharp peak even after the cancellation of the bulk and surface contributions. A 
peak associated with u = 2 A still shows up. As mentioned above, for q» > 0.025A/-uf, the 
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minimum of F results in a very broad maximum at u/2A m 0.8. We also note that for 
q\\ ^ 0.025A/t>^, while both the bulk and surface contributions are roughly proportional to 
q 2 , the net spectrum after the mutual cancellation of bulk and surface contribution effectively 
decreases as q\\ increases. 

For c/S, q\\c, and k±c <C 1 appropriate for layered superconductors, the total Raman 
intensity given in (|H]) can be reduced after some calculation to 



K J 2c V F 



2 



+ (2k ± c) 



(14) 



which is valid for qu <C u> < 1 and for all values of a. We note that in (14]), the function F 
given in is a very sensitive function of c/S and k±c and therefore we cannot approximate 
it by (u — l) 2 in the limit of c/8, k±c — > 0. We can see directly from (|TJ]) that the poles 
of b = ±1 are removed as a consequence of the cancellation between bulk and surface 
contributions. The only pole now is given by F = 0. One has a broad spectrum without 
any sharp peak for q\\ > 0.025A/-yp. In this region, one may verify that F oc E 2 + oc qi. As a 
result, the net Raman intensity I{oj) decreases as q\\ increases, roughly proportional to q^ 2 
(see Fig. ||). 

In addition to the surface contribution discussed above, there may be a surface plasmon in 
the region u > 2 A (see Fig. [|) but this only arises when we have different dielectric constants 
inside and outside the superlattice (e ^ e ) • In Fig. || we show the contribution of a surface 
plasmon to the Raman intensity. The parameters used are as in Fig. |3| but at a much higher 
momentum transfer q» = 0.2A/vp (which is probably the upper limit for Raman scattering 
in high-T c superconductors). In order to have a surface plasmon energy fairly close to 2A, 
the dielectric constants are taken to give a = —0.80. The latter value requires eo/e = 9, 
i.e., a layered superconductor with much lower dielectric constant compared to the overlay 
material. As shown in Fig. ||, this gives a Raman spectrum with lots of structure, with a 
broad surface plasmon peak at an energy above 2A. This region is the pair-breaking region, 
where there is strong BCS particle-hole damping of collective modes. The intensity of the 
surface plasmon is roughly proportional to q 2 . Once again, Fig. ||] shows the almost complete 
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cancellation between the bulk and surface contributions in the region of u < 2A. This is 
because at this relatively large value of gy, there is no well-defined solution of F = (see 
Fig. 0). 

One might worry that the bilayer out-of-phase phase mode shown in Fig. [1] might also 
be strongly modified due to the surface contributions at higher values of q». However, while 
the intensity (and the dispersion relation) of the out-of-phase phase mode is very sensitive 
to the pair tunneling strength Tj, it is not strongly dependent on the value of q\\. It is 
associated with out-of-phase oscillation of order parameters in a single bilayer, with no net 
charge fluctuation [p]]. Therefore, we have no reason to expect any strong modification of 
the out-of-phase phase mode at larger values of q\\. 

IV. CONCLUSIONS 

We have shown that the surface contribution plays a major role in determining the final 
Raman-spectrum line shape from semi-infinite superconducting superlattices. We find that, 
as in the case of normal superlattices discussed by Jain and Allen [|J, the proper inclusion 
of the surface contribution cancels "spurious" bulk contributions associated with Van Hove 
singularities (see Ref. 0) of the upper and lower limits of the bulk plasmon band of an 
infinite superconducting superlattice. We have also found that the bulk plasmon mode in 
the region u < 2A ceases to be well-defined when q» reaches a critical value (~ 0.025A/vp 
for the parameters we have used). In addition, the surface plasmon usually contributes to the 
Raman intensity in the region well above 2 A and only approaches 2A if the superconductor 
is overlayed by a transparent material with a dielectric constant eo much larger than that of 
the superlattice e. The surface-plasmon mode intensity increases as g| and thus one wants q\\ 
as large as possible if one wants to study it (see Fig. |6[). As mentioned above, for relatively 
large values of q\\ , the (negative) surface contribution tends to completely cancel out the bulk 
contribution in the entire frequency region below 2 A. Because the out-of-phase phase mode 
discussed in Ref. IjJ is not expected to be strongly affected by the surface contribution, at 
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relatively large values of q\\, this mode (see Fig. |]) should be the only structure remaining 
in the Raman spectrum below 2 A. 

The absolute intensity of the Raman spectra we discuss in this paper is somewhat below 
current experimental sensitivity. However we hope that the interesting predictions we make 
concerning the out-of-phase phase modes as well as surface plasmons will encourage future 
experimental efforts. 
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FIGURES 



FIG. 1. The bulk Raman intensity [11| given by (||) for a semi-infinite superconducting bilayer 
superlattice with interlayer Cooper-pair tunneling of amplitude Tj. We use Tj = 0.03g (x = 0.03) 
and gy = 5.0 x W~ 5 A/vf- The surface contribution is not included. 

FIG. 2. The dispersion relation of surface plasmons in a semi-infinite superconducting super- 
lattice (solutions of Q = 0) for a = —0.8 and —0.9. The shaded area between b = — 1 and u = 2 A 
lines denotes the bulk plasmon band for an infinite superconducting superlattice. The F = line 
represents a sharp "bulk plasmon" picked up by Raman scattering. For gy > 0.025A/fi? (dashed 
line), there is only an over-damped resonance (see Fig. ||). 



FIG. 3. Raman intensity given by (10) from both bulk (dotted line) and surface (dashed line) 
contributions for a semi-infinite superconducting superlattice. The same parameters are used as 
in Fig. |], with a = 0.82 (corresponding to eo = 1 and e = 10). Apart from the absence of the 
out-of-phase phase mode (only expected when the unit cell has two layers), the bulk contribution 
is very similar to that given in Fig. [I]. 

FIG. 4. Raman intensities as in Fig. |3| except for a much larger value q» = O.lA/vp. This plot 
clearly shows the almost complete mutual cancellation of the bulk and surface contributions at 
such large values of q\\. 

FIG. 5. The total Raman intensity in the region u < 2A (including both bulk and surface 
contributions) for various values of the momentum transfer gn. The peak is associated with the 



zero or minimum of F in (14) 



FIG. 6. The Raman intensity as in Fig. ^, with qn = 0.2A/f_p and a = —0.8. 
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